Midgap States in Antiferromagnetic Heisenberg Chains with A Staggered Field 
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We study low-energy excitations in antiferromagnetic Heisenberg chains with a staggered field 
which splits the spectrum into a longitudinal and a transverse branch. Bound states are found 
to exist inside the field induced gap in both branches. They originate from the edge efi'ects and 
are inherent to spin-chain materials. The sine-Gordon scaling h?s^^\\ogha\^^^ {hg-. the staggered 
field) provides an accurate description for the gap and midgap energies in the transverse branch for 
S = 1/2 and the midgap energies in both branches for S — 3/2 over a wide range of magnetic field; 
however, it can fit other low-energy excitations only at much lower field. Moreover, the integer- 
spin S — 1 chain displays scaling behavior that does not fit this scaling law. These results reveal 
intriguing features of magnetic excitations in spin-chain materials that deserve further investigation. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.Mg, 75.40.Gb 



Quantum spin chains have been a source of fascinat- 
ing developments in physics and manyrelated fields since 
the early days of quantum theory They serve 

as model systems for revealing fundamental physics in 
many materials. Progress in experimental techniques has 
led to the synthesis and characterization of many quasi- 
one-dimensional materials in the last two decades. A 
very active field of theoretical and experimental research 
has emerged for the study of novel magnetic properties 
in low-dimensional spin systems where quantum fluctua- 
tions play a crucial role 0, [|| . 

Intriguing low-energy phenomena like the spin 
gap induced by magnetic field have been observed 
in many recently synthesized materials such as 
Cu(C6D5COO)23D20 0, iK Yb4As3 Til and 
CuCl9 2 fdimethvlsulfoxide) jl7| . The Dzyaloshinskii- 
Moriya (DM) interaction is intrinsic to these materials 
due to an alternating structure of molecular or spin- 
orbital interactions ^3, ^3- Although it is one order 
of magnitude smaller than the standard Heisenberg ex- 
change interaction in these materials, significant effects 
of the DM interaction have been unveiled recently when 
the materials are subjected to external magnetic fields 
0, 1^ Il4j . An energy gap proportional to some power 
of the magnetic field opens up. These systems can be 
generally modeled by the following Hamiltonian, 



H 



JS, • S,+i - (-)'D • S, X S,+i 



(1) 



where the three terms in the summation are the antiferro- 
magnetic Heisenberg, DM and Zeeman splitting interac- 
tions, respectively. While the exchange coupling constant 
J and the uniform (staggered) g" (g^) can be determined 
from the neutron scattering and electron spin resonance 



experiments, the D-vector can be obtained only through 
theoretical analysis of the gap values in comparison with 
experiment |3. ImIi5I |. 

The DM term can be eliminated by performing a 
spin rotation about the D vector by an angle a = ±i 
tan~^(Z?/J) on the alternating sites. After neglecting 
anisotropic terms in the limit of I? <C J, which are only 
needed to account for the field dependence of the gap 
in different crystallographic directions T^, the relevant 
physics is captured by the effective Hamiltonian 



H = J^[s, • s.+i - hjf - i-yhj^ 



(2) 



where hu and hs are the effective dimensionless (scaled 
by gfis/J, and J is taken as the energy unit hereafter) 
uniform and staggered field, respectively. This Hamil- 
tonian has been mapped onto the quantum sine-Gordon 
model using the bosonization technique 4] and the afore- 
mentioned novel magnetic properties are described in 
terms of soliton, antisoliton, and breather excitations 

Despite intensive investigations, quasi-one-dimcnsional 
spin systems with DM interaction continue to display sur- 
prising new physics. In this paper we show that the ex- 
citation spectra are split into a longitudinal and a trans- 
verse branch, while edge effects result in midgap states 
in both branches for open chains with both integer and 
half-integer spins in the thermodynamic limit. These 
midgap states are inherent to real spin-chain materials 
which always have chain ends. The existing field the- 
ory is established in the context of the quantum sine- 
Gordon model for periodic spin- 1/2 chains 0i|l3- It is 
very important to clarify how soliton, antisoliton, and 
breather modes in the field theory are related to the gap 
and midgap states in different branches. Recent elec- 
tron spin resonance experiment |l8l | on a spin- 1/2 chain 



2 



compound [Pyrimidine-Cu(N03)2(H20 )2]„(CuPM) re- 
vealed that some excitation modes, including one of the 
most intensive modes, cannot be fuUydescribed in terms 
of the quantum sine-Gordon model 4]. This raises fun- 
damental issues regarding the nature of these magnetic 
excitations. Moreover, current discussions focus almost 
exclusively on the S = 1/2 case; it is of great interest to 
examine the low-energy excitations of higher-spin chains 
[l^ Efil | and their possible connections to the spin- 1/2 
chain and the field theoretical description. 

We use the densit y m atrix renormalization group 
(DMRG) method Hlli^l to study directly the Hamil- 
tonian defined by Eq. (j^J for both integer (S'=l) and 
half-integer (5'=l/2 and 3/2) spin chains. We find the 
same qualitative behavior for the gap and midgap states 
over a wide range of and, consequently, focus below 
on the effect of hs at hu=0. Unless explicitly stated oth- 
erwise, we use open boundaries in all calculations. The 
results are obtained by keeping 400 states, but exam- 
ined at low fields with 800 states. The truncation errors 
are on the order of 10^^ for S = 3/2 and much smaller 
for 5*= 1/2 and 1. The gap and midgap energies in the 
thermodynamic limit were obtained by extrapolating nu- 
merical results for up to 200 sites to the long-chain limit. 
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FIG. 1: The staggered field dependence of the transverse 
and longitudinal gaps (At: solid circles; A^: solid triangles) 
and the midgap-state energies (A™: open circles; A™: open 
squares) for the spin-1/2 chain. The lines are fittings to the 
hi^^\^oghs\^^^ scaling behavior anticipated from the confor- 
mal field theory The slopes of the plots which determine, 
up to log corrections, the scaling exponents of the leading or- 
der term are shown in the inset with the same axis labels, in 
the logarithmic scale. Data show a common scaling exponent 
of 7 = 2/3 when the field is not too high. 

Figure ^ shows the low-lying energies relative to the 
non-degenerate ground state for the spin-1/2 chain. At 
hu — the spectrum splits into a transverse branch 
with \S^^^\ = 1 and a longitudinal branch with S^g^=Q, 
with gaps At and A^, respectively, determined by scal- 
ing analyses that confirm the continuous nature of each 
spectrum in the large-L limit. Midgap states appear in 



both branches. The gaps and the midgap-state energies 
A™ and A™ are obtained from the extrapolations of the 
DMRG data as shown in Fig. 2. 




0.01 0.02 0.03 0.04 0.05 
1/L 

FIG. 2: The excitation energies versus 1/L for several lowest 
states of the spin-1/2 chain labeled by 3^*°* (upper panel) 
and Sf'"' {\Slot\) (lower panel) [a = A, B, C]. They all scale to 
zero as L — > oo at /is =0. At finite hs (hs — 0.05 shown), they 
split into a longitudinal (lines) and a transverse (symbols) 
branch with distinct gap and midgap states (^^(l) — > A™, 
SUl),Sh{l) ^ At, S\{0),Sh{0) ^ A^, S^(0) ^ A^, see 
the text for details). 

The conformal field theory for the spin-1/2 chain un- 
der a staggered magnetic field predicted 4] the scaling 
behavior hl^^llog hs\^^^ for the induced gap. We apply 
this scaling form to our numerical results as shown in 
Fig. 1. It is interesting to note that the field theory 
not only provides a good description in the low-field re- 
gion for the transverse gap, a periodic chain feature, as 
expected, but also fits well the transverse midgap, a gen- 
uine open boundary edge effect (see below). The scaling 
exponents at the low field limit share the common value 
of 2/3 predicted by the field theory. However, the fitting 
range is different for the two branches: the deviation from 
the predicted scaling behavior starts at hg « 0.03 in the 
longitudinal branch while a satisfactory fitting persists 
beyond hs=0.1 for the transverse branch. Moreover, the 
midgap state appears in the transverse branch as soon as 
hs is applied, but it does not become distinguishable in 
the longitudinal branch until hg =0.1 which corresponds 
to a field of 14 Tesla in copper benzoate 1^. It suggests 
that the midgap state in the transverse branch should 
be observable at low fields, while much higher fields are 
required for its identification in the longitudinal branch. 
Since the midgap states are due to the chain-end effect. 
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their intensity is proportional to the impurity /defect con- 
centration. Samples with high level of nonmagnatic im- 
purity doping may be needed for their detection. 

To examine the origin of the midgap states, we consider 
a bond impurity model |2^. A coupHng J' is added be- 
tween the two end spins and varied from (open bound- 
ary) to J (periodic boundary). The midgap states de- 
velop in each branch when J' deviates from J, indicating 
that these states appear as a consequence of the edge 
effect. The same effect is also observed in all higher- 
spin chains discussed below. Alternatively, we show in 
Fig. 2 for the spin-1/2 case that the midgap states can 
be traced back to the three lowest excitation states with 
Slot = 1 at /is=0, labeled as S\, and S'^, respec- 
tively. A comparison of the /is = and 0.05 data shows 
that state S\ splits into midgap states S\{1) and S'^(O), 
state splits into ^^^(l) that becomes the bottom of the 
transverse continuum and S'^{Q) that is a midgap state 
degenerate with S\{Q), and state splits into S}j{\) 
degenerate with S^{1) and 5*^(0) which is the bottom of 
the longitudinal continuum. The lowest excitation state 
S\ with Stot = moves higher in energy and enters the 
longitudinal continuum at /i<; > 0, which mixes two states 
S\{0) and 5^(0) at about L = 36. 

We now turn to the spin-1 case where the spectrum 
is gapful at /is = 0. In the thermodynamic limit, the 
ground state of the open chain is four-fold degenerate 
owing to the topological edge effect as interpreted in 
the valence-bond-sohd (VBS) picture The stag- 

gered field splits the four-fold degenerate ground state 
into three mixed states: a S^^^ = ground state, a dou- 
bly degenerate \S^g^\ = 1 transverse midgap state and a 
S'fot = longitudinal state. The dependence of the gap 
and midgap-state energies on hs is shown in Fig. |31 In 
addition to the midgap states originating from a mixture 
of the \S^^^\ = 0, 1 states in the ground-state manifold, 
there are also midgap states with S^^^ — and \S^g^\ ~ 2 
below but close to the bottom of the longitudinal contin- 
uum. We have tried fitting the results to the field theo- 
retical formula, but found that the /is^'^l \oghs\^^^ scahng 
is unsuitable for S* = 1 in any field range. While this is 
not surprising for the gaps at low field since they have 
finite values at hs = 0, the fact that it does not apply 
to the midgap energies that do scale to zero at /i^ = 
strongly suggests different intrinsic behavior for integer 
and half- integer spin cases |3|. 

Figure 01 shows the scaling behavior of the gap and 
midgap-state energies of the spin-3/2 chain. The low- 
field sine- Gordon scaling derived for the spin-1/2 chain 
also fits these results. In particular, the two lowest 
midgap states are well described berond /is = 0.1. It 
highlights common central charge j25| and topological 
features intrinsic to all half-integer-spin chains. How- 
ever, subtle yet important differences exist between the 
S" = 1/2 and S' = 3/2 cases. Unlike the spin-1/2 case 
where bound midgap state develops at low field explic- 
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FIG. 3: The staggered field dependence of the transverse and 
longitudinal gaps (At: solid circles; Al: solid squares) and 
the midgaps (A™: open circles, \Stot\ = 1; A™: open squares, 
Stot ~ 0; down-triangles: Stot ~ 0, diamonds: Sfot = 0, and 
up-triangles: \Stot\ = 2) in the spin-1 chain. The lines are 
guide to the eyes only. The midgap states falling off from the 
longitudinal continuum are not connected by lines for clarity. 



itly in the transverse branch only, midgap states appear 
in both branches at very low fields for S = 3/2. More- 
over, in the spin-3/2 case the /is^^ | log /is 1^^^ scaling fits 
both gaps only in the low field limit (/ig < 0.02), while 
the fitting works very well for the two lower midgap states 
in both branches over a much larger field range. Mean- 
while, similar to the S=l case, additional midgap states 
fall off from the bottom of the longitudinal continuum 
(not shown in Fig. 4 for clarity since a few of them are 
very close to the bottom at low fields). These similar- 
ities and distinctions between different half-integer spin 
chains deserve further studies. 

Finally, we remark on the nature of the midgap states 
that share some common topological origin for all spin 
magnitudes. Although the gap at /is = is considered 
a fundamental criterion to distinguish the integer and 
half-integer spin cases 0, Q j it has been shown that an 
extra Berry phase contribution reflecting the topologi- 
cal feature leads to edge states in all open spin chains 
with S" > 1/2 Most recently, it was also shown 2^ 
that a topological string order originally derived |2^ . 29j 
for integer-spin chains exists in half-integer-spin chains 
as well. Nonetheless, the topological edge effect has been 
experimentally identified so far only in integer-spin chain 
materials where the excitation spectrum is gapful and 
topological excitations unstable against external distur- 
bances (such as doping), can be measured [i^. For the 
S = 1 case, this effect was observed in Yi_a;BaCaa;Ni05 
[30| and Y2BaNii_2,Mg3;05 [U- The present results un- 
veil the topological edge effects driven by the staggered 
magnetic field. The midgap states split off from the de- 
generate ground state of S" = 1 chain (and all other 
integer-spin chains) represent a secondary (field induced) 
topological edge effect since bound midgap states are al- 
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ready observable at /i^ = 0. Meanwhile, additional bound 
states fall inside the field induced gap from the contin- 
uum spectrum of all spin chains with S > 1 at > 0. 
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FIG. 4: The staggered field dependence of the transverse and 
longitudinal gaps (At: solid circles; Al: solid squares) and 
the corresponding lower midgap-state energies (A™: open cir- 
cles; A™: open triangles) in the spin-3/2 chain. The lines are 
fittings to ht^^ I log hs \ for these low-energy excitations. 

It should be noted that the staggered field hg arises 
from a nonzero uniform field hu in real materials. The 
latter does not change the qualitative scaling behavior 
when the field is not too strong 0, We have ex- 

amined the gap and midgap states at nonzero and 
found the same qualitative results at realistic hg/hu ra- 
tios although the classification of the longitudinal and 
transverse branches no longer holds strictly because of 
the breaking of the axial symmetry by 

In summary, we have shown that bound midgap 
states generally exist in open boundary antiferromagnetic 
Heisenberg chains with a staggered magnetic field. They 
should appear in quasi-one-dimensional materials with 
the Dzyoloshinskii-Moriya interaction. Some of the gap 
and midgap energies for the half-integer spin chains fit 
well to a scaling function derived from the quantum sine- 
Gordon model. However, it is revealed that (i) other low- 
energy excitations of the half-integer spin chains do not 
fit equally well and (ii) the scaling behavior of the integer 
spin chain is qualitatively different. Further experimental 
and theoretical investigations are called for to fully iden- 
tify low energy excitations, especially the midgap states 
as a general phenomenon in open spin chains. 
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